On Time Delay Margin Estimation for Adaptive Control and 

Optimal Control Modification 

Nhan Nguyen* 

NASA Ames Research Center, Moffett Field, CA 94035 


Abstract 

This paper presents methods for estimating time delay margin for adaptive control of input delay systems 
with almost linear structured uncertainty. The bounded linear stability analysis method seeks to represent an 
adaptive law by a locally bounded linear approximation within a small time window. The time delay margin 
of this input delay system represents a local stability measure and is computed analytically by three meth- 
ods: Pade approximation, Lyapunov-Krasovskii method, and the matrix measure method. These methods are 
applied to the standard model-reference adaptive control, cr-modification adaptive law, and optimal control 
modification adaptive law. The windowing analysis results in non-unique estimates of the time delay margin 
since it is dependent on the length of a time window and parameters which vary from one time window to the 
next. The optimal control modification adaptive law overcomes this limitation in that, as the adaptive gain 
tends to infinity and if the matched uncertainty is iinear, then the ciosed-loop input deiay system tends to a 
LTI system. A iower bound of the time delay margin of this system can then be estimated uniquely without the 
need for the windowing analysis. Simulation results demonstrates the feasibility of the bounded linear stability 
method for time delay margin estimation. 


1 Introduction 

Input delay systems are generally non-minimum phase. For linear input delay systems, feedback gain must be kept 
to a reasonable value to maintain stability. Input delay influences stability of adaptive control in a similar manner. 
Adaptive gain is used to control the rate of adaptation in adaptive control. For model-reference adaptive control, 
it is well-known that as the adaptive gain increases, the closed-loop system loses robustness, thereby rendering it 
susceptible to instability in the presence of unmodeled dynamics and or input time delay. Thus, to maintain stability 
of an input delay adaptive system, the adaptive gain must be carefully selected. For a given value of the adaptive gain, 
there exists a corresponding value of input time delay for which the adaptive system is on the verge of instability. This 
is known as a time delay margin. To maintain stability, the adaptive gain of the system must be kept below the value 
that corresponds to the time delay margin of the system. 

Global stability analysis for input delay adaptive systems is a challenging problem. Lyapunov-Krasovskii method 
or Lyapunov-Razumikhin method are much more difficult to apply to an adaptive system. Even for a simple scalar LTI 
system, both the Lyapunov-Krasovskii method or Lyapunov-Razumikhin method can result in conservative estimates 
of the time delay margin [1], So, even if a global stability analysis for input delay adaptive control is available, the 
conservatism in the estimation could render it impractical. The lack of available analytical methods for computing the 
time delay margin of adaptive control is a hurdle for certification of adaptive control [2], 

While global stability analysis is challenging, several studies have recently been done to address local stability of 
input delay adaptive systems. One such method applies a Pade approximation to transform an input delay system into 
a delay-free higher-order system [3]. The transformed system is then analyzed using the standard Lyapunov method 
to estimate the time delay margin. However, this approach yields a highly conservative time delay margin even for a 
simple scalar adaptive system [3] thereby rendering the approach rather less practical. Another method, called bounded 
linear stability analysis, attempts to analyze the stability margins of an adaptive system in a local context [4, 5]. The 
method approximates an adaptive system as a series of bounded linear systems inside time windows. The windowing 
analysis allows the bounded linear systems to be analyzed using linear analytical tools. The method has been shown 
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to provide a less conservative estimate of the time delay margin. Subsequently, a similar method has been developed 
using the windowing approach to estimate the stability margin of an adaptive system [6]. LMI methods also has been 
used to analyze stability of adaptive control [7]. The Lyapunov-Razumikhin method has been used to estimate the time 
delay margin for a simple scalar adaptive system [8]. The method requires optimization of the candidate Lyapunov 
function in order to reduce the conservatism in the estimated time delay margin. 

This paper extends the bounded linear stability analysis method for analyzing input delay adaptive control with 
almost linear structured uncertainty. Stability of an input delay adaptive system is analyzed by three methods: Pade 
approximation, Lyapunov-Krasovskii method, and the matrix measure method, to estimate the local time delay margin 
of a bounded linear system inside each time window. Three different adaptive laws are used: the standard MRAC, 
cr-modification adaptive law [9], and optimal control modification adaptive law [10]. Asymptotic analysis of the time 
delay margin as the adaptive gain tends to infinity is performed to study to effect of large adaptive gain on the time 
delay margin. Simulations are studied to demonstrate the feasibility of the time delay margin estimation. 


2 Input Delay Adaptive Systems and Bounded Linearity Stability Analysis 


Given an input delay nonlinear plant 


x(t) =Ax(t) +B u(t — td) + ©* T <l>(x(f)) 


( 1 ) 


where x(t) : [0,<=°) — > R" is a state vector, u(t) : [0,°°) — > R p is a control vector, A G R" x " and B G W xp are known 
such that the pair (A . B) is controllable, 0* G R mxp is an unknown constant weight matrix that represents a parametric 
uncertainty, <t> (. x (t)) : R" — > R m is a vector of known functions, and tj is an input time delay. 

The structure of the uncertainty is assumed to be linearly dominant. That is 

<3>(x(f)) = x(t) + 8 (x) (2) 


where ||5 (x)|| <C ||x(f)|| is small. 

The input delay tj could also be viewed as the time delay margin for robustness against unmodeled dynamics of 
the delay-free system 


x(t)=A\(t)+B m (/) + ©* <3>(x(f)) 


(3) 


The reference model is specified as 

Xm (t) = A m x m (?) T B m r(t ) 


(4) 


where A m G R” x ” is Hurwitz and known, B m G R' ,xp is also known, and r (t) : [0,°°) — » R p G Jzfoo is a bounded 
command vector. 

Defining the tracking error as e ( t ) = x m ( t ) — x(f), then the controller u ( t ) is specified by 


u ( t ) = -K x x ( t ) + K r r (t) - u a d (t) 


(5) 


where K x G R px ” and K, G R pxp are known nominal gain matrices, and u a d (t) G R p is an adaptive signal given by 


Had (t) = © T (t) & (x) 

where 0 G R mxp is an estimate of ©*. 

Assuming that the model matching conditions can be satisfied, then 

(6) 

A m =A — BK X 

(7) 

B m = BK r 

The standard model-reference adaptive control law is 

(8) 

©(f) = -rd>(x(f))e T (t)PB 

where T = T t G W ixm > 0 and P = P G R' !X,i > 0 solves the Lyapunov equation 

(9) 

PA m +Aj n P = — Q 

(10) 
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where Q = Q r G R” x " > 0. 

The error equation corresponding to the input delay system (3) can be derived by substituting the time-delay 
version of the controller from Eq. (5), thus resulting in 

e (t) = Ae (t) - BK x e (t - t d ) + Bu ad (t - t d ) - B &* T (x (t)) - BK X [x m (, t ) - x m (t - t d )\ + BK, [; r (t) -r(t-t d )\ (11) 

To analyzing this system, the bounded linear stability analysis method has been proposed to approximate the 
adaptive system as a bounded linear system within a given time window [4, 5]. The windowing analysis then permits 
the use of linear tools to analyze stability of the approximated bounded linear system inside a given time window. 
Theorem 1: The adaptive law (9) is bounded locally by a linear approximation as 

© T (f)4>(x(f)) = -yB J Pe(t) (12) 

where y is a constant defined locally and retrospectively as 

7= Jr [ $ T (x(T))r<E(x(T))c/T (13) 

T ° Jtt-To 

for t £ [tj — where to = 0, ?, = f,-_ ] + Tq and i = 1,2 ,n 

Proof: Choose a Lyapunov candidate function 

V (t) = e T (t)Pe (t) + trace [<g) T (r)r 1 G)(r)l (14) 

The error equation of the delay-free system is 

i (f) = A m e (f) + B& J (t) T> (x (f)) (15) 

where @ (t) = 0 (t) — &* . 

Denoting V g (t) as the Lyaponov candidate function to be evaluated globally using the adaptive law (9) as follows: 

V g (t) = -e J {t) Qe(t) + 2e J (t)PB& J (f)4> (x(f)) - 2trace © T (f)4>(x(f))e T (t)PB = -e T (f) Qe (t) < 0 (16) 

Denoting V/ ( 1 ) as the Lyapunov candidate function to be evaluated locally with a time window using the locally 
bounded linear approximation (12) yields 


Vi ( t ) = — e 1 (t)Pe(t) +2e T (t)PB@ T (t)<P(x(t)) — 2trace 


= V g (t) + 2e r ( t)PB@ T (t)0(x(t)) — 2trace 


T>(x(f))-2trace |© T (^r -1 ©^) 

— 2trace|© T (t)<P(x(t)) <t> T (x(t))r<P(x(t)) <t > 1 (x(t)) @(f)| 

y g (f)+2e T (f)PB© T (f)T>(x(f)){l-7 4> T (x(f))r4>(x(f)) | (17) 


for t £ [tj — To,tj), where to = 0, tj = tj _ i + To and i = 1 , 2 , . . . , n 
Consider the integral form of Eq. (17) 


rU ftj rtt 

/ Vj (t) dt= V g (t) dt+ 2e T 

Itj-To Jtj—To Jtj—To 


(t)PB& r (f)4>(x(f)) |l — y <f> T (x(f))r<I>(x(f)) l }dt (18) 


Then for V g ( t ) and V/ (t) to be equal in a time window 


Vj\t)dt= / V g (t)dt 


which implies 


(t)PB© T (f)T>(x(f))|l -7 T> T (x(f))r<I>(x(f)) l }dt = 0 
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It is important to note that this is a definite integral equation for which a valid solution can include a constant 
solution of y. Such a solution is called a “weak-form” or integral-form solution which is valid only over a finite 
time interval. In contrast, the “strong-form” solution is a global solution that satisfies for all time. In the windowing 
analysis, the weak-form solution is used. 

The mean value theorem for integration states that 

rb rb 

/ F(t)G(t)dt = F(c) / G(t)dt (21) 

J a J a 

where c £ [a,b] and G(?) > 0. 

Let t £ [; tj — T 0 ,ti), then applying the mean value theorem for integration to Eq. (20) yields 


ru r i - 1 r 

/ 2e J (t)PB& J (?)<t>(x(?)) <t> T (x(t))r<P(x(t)) <t> T (x(t))r<t>(x(t))-y 
>ti-T 0 1 J L 


dt = 


+ 2 e T (1) PB& J (?) <t> (x (?) ) | <t> T (x (?) ) DI> (x (?) ) 

Hence, (13) is thus obtained. Then it follows that 




<P r (x(t))r<P(x(t))dt~YT Q 


= 0 (22) 


V, (i t; ) - V, ( ?; - To) = V g (, ti ) - Vg (?; - To) < 0 (23) 

Thus, the local Lyapunov candidate function V) (?) is a piecewise approximation of the global Lyapunov candidate 
function V s (?) where their values are equal at the beginning and end points of a time window. 

Using the bounded linear approximation of the adaptive law (9), one gets a piecewise locally bounded linear 
approximation of the standard MRAC adaptive law (9) 

*ad(t) = © T (?)<L(*(?)) + 0 T (?)4> (*(?)) « -jB T Pe(t) + & T (?)<*>(*(?)) (24) 

for ? £ [tj — To,?,-), where ?o = 0, ?; = tj _ i + To, and i = 1 , 2, . . . ,n —> °°. 

The second term in the right hand side can be locally approximated by a first-order Taylor’s series as 


© T (?)<i>(*(?)) = 0T£ 

T=i 


d^{x(ti)) . , , T XT' <? 2 <L(x(?;)) . , . . . , 

X j (tj) +0; 52 52 3 ,, 3 .... X j fa) \- Xk M Xk ( f ')] 


dxj ^ 


j^ltl dx J dxk 


•©; £ fAI> ^ j,/ii ;.v ; i?)-.v ; i? ) i • ... 


j= i 


= ©2 [x m (?) - e (?)] + ©2 ***,- Xi [x m (?) - e (?) - Xj] + ... (25) 


where ^$'2) = <5<& (x (tj)) /dxj and (&'■ ' x^j = E^ =1 x k (tj) cf<l> (x (tj)) /dxjdx k , j=l,...,n. 

Thus, the locally bounded linear approximation of the error equation and the standard MRAC adaptive law are 
expressed as 


e (?) = Ae (?) - BK x e (t-t d )+ Bu ad (t - t d ) - B0* T <I>; - B0* T <I>'. [x m (?) - e (?) - Xj) - BK X [x m (?) - x m (t-t d )\ 

+ BK r [r(t)-r(t-t d )\ (26) 

u ad (t) = ~yB T Pe (?) + ©2 \x m (?) - e (?)] (27) 

for ? £ [tj — To,?;), where ?o = 0, tj = ?;_i + To and i= 1,2, . . . ,n — > °°. 

Equations (26) and (27) show that the stability of the locally bounded linear approximation depends on several 
factors: 

• The initial condition x (0) 

• The structure of the matched uncertainty <1> (x (?)) 

• The parametric uncertainty 0* 
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• The plant and reference model matrices A , B, A m , and B m 

• The adaptive gain parameter y which includes the adaptive gain T as well as the square of the amplitude of 

<W0)- 

Thus, it can be seen that this bounded linear approximation does capture the complex nature of stability of a nonlinear 
adaptive control system, at least in a local sense. 

In a special case when <t>(x(f)) = x(t), then the bounded linear approximation of the error equation and the 
standard MRAC adaptive law become 

e (t) = Ae(t) — BK x e ( t-t d ) + Bu ad ( t-t d )~ B&* T [x m (t) -e(t)\- BK X [x m (t) - x m ( t - t d )\ + BK r [r (f ) - r {t - t d )\ 

(28) 

u (l d it) = - JB T Pe ( t ) + [x m (t) - e (f)] (29) 

3 Time Delay Margin Estimation of LTI Systems 

Consider an input delay closed-loop LTI system 

x it ) = Ax it ) — BKx it — t d ) (30) 

where x(t) : [0,°°) — > K" and A iA — BK) £ C”, i.e., A — BK is Hurwitz. 

The time delay margin is defined by the following characteristic equation 

det ijcol- A +BKe- jmd ) =0 (31) 

For simple systems, analytical solutions of t d can be computed, but in general such solutions are not easily obtained. 
We present three methods for estimating the time delay margin. 


3.1 Pade Approximation 

The Laplace transform of the input delay LTI system is 

(s) = AX (s) - BKX (s) e~‘ dS 

Consider the following first-order Pade approximation 

e -t d s _ 2 ~tdS 
2 + t d s 

Then the approximate input delay system becomes 

(2 + t d s) sX (s) = (2 + t d s) AX (s) — (2 — t d s) BKX is) 
In the time domain, this is expressed as 

tdxit) = (— 2I + t d A + t d BK)xit) +2 (A — BK) x(f) 


The time delay margin is then found by 


det 


(0 2 I + jco I+A + BK H — (A — BK) 


td 


td 


= 0 


Alternatively, the time delay margin can also be obtained as 


det 


jot! 


- 1 


-UA-BK) ja+ll-A-BK 


= 0 


(32) 

(33) 

(34) 

(35) 

(36) 

(37) 
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Example: Given 


A = 


0 1 

-1 1 


BK = 


0 0 
0 2 


The time delay margin using the first-order Pade approximation is estimated to be t d = 0.528 sec. 
The exact results can be determined from 


det 


j(0 - 1 

1 ja> — ] + 2 (cos m* d — j sin cot d ) 


= — co 2 + 2cosincot d + 1 — jco(l — 2coscot d ) = 0 


V3 + V7 

co = Z =2.189 rad/sec 


td = 


2 n 


= 0.478 sec 


3(^3 + V7) 

So the time delay margin by Pade approximation is non-conservative. 


3.2 Lyapunov- Krasovskii Method 

Stability of time-delay differential equations based on the Lyapunov-Krasovskii method using Lyapunov-Krasovskii 
functionals have been studied exhaustively [11, 12]. The negative-definiteness of the time derivative of a Lyapunov- 
Krasovskii functional results in a linear matrix inequality that can be solved for a time delay margin. The solution is 
generally non-unique since it depends on the choice of a Lyapunov-Krasovskii functional. As a result, the time delay 
margin obtained by the Lyapunov-Krasovskii method can be conservative. 

Theorem 2: Lor the closed-loop input delay system (30), the system is asymptotically stable if the following linear 
matrix inequality is satisfied: 

(A - BK) J P + P (A - BK) + (a + j3 ) t d PBKK T B T P + ^A T A + jK t B t BK <0 (38) 

The time delay margin is the largest value that renders the LMI feasible. 

Proof: Lor the input delay LTI system (30), we write 


Then 


/ x(x)dx = x(t) — x(t — t d ) 

Jt-t d 

x (t — t d ) =x(t) — f Ax(x)dx+ f BKx(x-t d )dx 

J t-t d J t-t d 


The input delay system now becomes 


x(t) = (A-BK)x(t)+BK f Ax(x)dx — BK f BKx{x-t d )dx 

J t-u J t-t d 


t-t d J*-t d 

Consider the following Lyapunov-Krasovskii functional [11] 

rt pt pt pt 


where P = P T >0. 

Evaluating V ( t ) yields 


(39) 


(40) 


(41) 


V(t)=x T (t)Px{t)+ f [ x T (e)Qx(e)d6dx+ [ f x T {6)Sx{e)dedx>0 (42) 


V (t) =x T (t)Px(t) +x T (t)Px(t) + f x T (t) Qx(t)dx — f x T {x)Qx{x)dx+ f x T (t — t d )Sx(t — t d )dx 

dt~t d J t~t d J t~ t d 

- [ x T (x-t d )Sx(x-t d )dx (43) 

Jt-t d 
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This becomes 


V (?) = x T (?) (A — BK) J P + P (A — BK ) x (?) + 2x t (?) PBK / Ax (r) c/t - 2x T (?) PBK / BAx (t - t d ) dx 


>t-t d 


>t-t d 


+ t d x T (?) Qx (?) - f x T (x)Qx(x)dx+t d x T (t)Sx(t)~ f x T (x-t d )Sx(x-t d ) dx (44) 

Jt-t d Jt-t d 


From the following quadratic expressions 


at d x T (t)PBKK T B T Px(t)-2x J (t)PBK [ Ax(x)dx+ — ([ x T (x)A J dx]([ Ax(x)dx) > 0 (45) 

Jt-t d at d \Jt-t d ) \Jt-t d J 


pt d x T (; ? ) PBKK t B t Px (?) + 2x t (?) PBK f BKx {x - t d ) dx 

Jt-td 

+ W,(f xT ( t_ ^)^ TbTJt ) (/ BKx(x-t d )dx^ > 0 (46) 


we obtain 


2x J (t)PBK[ Ax(x)dx< at d x T (t)PBKK T B T Px(t) + — ( [ x T (x)A T dx]( f Ax(x)dx ' 

Jt-t d at d \Jt-t d J \Jt-t d J 

1 !' 

< at d x J (t)PBKK r B T Px(t) + - x T ( x)A J Ax(x)dx (47) 

a Jt-td 


- 2x t (?) PBK f BKx (x - t d ) dx < p t d x T (?) PBKK J B J Px (?) 

Jt-td 

+ (f xT ^ t d)K T B T dx^j ^ BKx(x-t d )dx 

1 

<pt d x T (t)PBKK T B T Px{t) + - / x T (T-? d )/s: T B T BKx(T-?,,)t/T (48) 

P 


for some a > 0 and p > 0. 

Choose Q=A J A/a and 5 = K T B T BK/p. Then 


V (? ) < x 1 (t)\(A-BK)' P + P (A- BK)\x (?) + at d x ' (?) PBKK J B J Px (?) + )3 t d x ' (?) PB/Of 1 B 1 Px (?) 


-t d t, 


+ |x T (?) A t Ax (?) + ^x T (?) K t B t BKx (?) (49) 


For stability, V (?) < 0. Thus we obtain the LMI (38). 

Example: Given 

’ 0 0 ' 

0 2 

Choose a = 1 and p = 1 . Let 

{A-BK) r P + P{A-BK) = -0 


A = 


0 1 

-1 1 


b/s: = 


and choose 2 = 7 where I is the identity matrix. Then the time delay margin is estimated to be t d = 0.068 sec. 

t d can be maximized by a suitable selection of a and p . Figure 1 is a plot of t d as a function of a and p . The 
maximum value of t d is 0.072 sec corresponding to a = 0.51 and p = 0.89. As expected, the Lyapunov-Krasovskii 
method produces a very conservative estimation of the time delay margin. 


7 




2 2 

a 

P 

Fig. 1 - Time Delay Margin by Lyapunov-Krasovskii Method 


3.3 Matrix Measure Method 


The matrix measure method has been developed recently and affords a simple way to estimate the time delay margin 
and the effective phase margin for the MIMO LTI system [1], The matrix measure fj. is defined as an eigenvalue of a 
symmetric part of a complex matrix [11] such that 


Hi (C) = A t 


/ C + C * 
V 2 


where C G C is a complex matrix and C* is its complex conjugate transpose, then /i has the following properties 

fli (C) G M 

||7+eC7||-l 
e 

1-p-eCll 
£ 

H (C) < ReA, (C) < ju (C) 

ImA (C) < Ji(-jC) 

Theorem 3: The input delay LTI system (30) is asymptotically stable if the following inequalities hold 


H(C) = 

max A j 

1 <i<n ' 

rc+c*\ 

| = lim 
1 e-+0 

V 2 ) 

H(C) = 

min A; 

\<i<n ' 

(c+c*\ 

| = lim 
1 e^O 

\ 2 ; 


1 

td < — cos 
CO 


H(A)+H{jBK) 


\\BK\\ 

co<Ji(-jA) + \\BK\\ 


where II. II = 


is the Jz? 2 -norm. 


(50) 

(51) 

(52) 

(53) 

(54) 

(55) 

(56) 

(57) 


Proof: The real parts of the system poles are bounded from above by 
a = ReA,- (A - BKe^ ja, ' d ) <Ji(A) + Jl (- BKe -> md ) <Ji(A)+Ji ( -BK ) |cos c 0 t d \ + Jl (jBK) |sin cot d \ (58) 

Let 0 < cold < then the input delay system is stable if <7 < 0 which implies 

A H A ) < —ji (—BK) cos cotd — H (jBK) sin <ot d = H (BK) cos cot d — ft (jBK) sin cot d 

< Jl (BK) cos cot d — Ji (jBK) sin cot d (59) 
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Upon some algebra, this can be expressed as 

[p 2 ( BK ) + p 2 ( jBK )] cos 2 cot d -2 p(A)p ( BK ) cos c Ot d +p 2 (A) -p 2 ( jBK ) > 0 (60) 

The solution yields a bound on time delay margin t d as 


t d < — cos 
co 


1 p(A)p(BK) + p(jBK) \Jp 2 (BK) +p 2 (jBK) - p 2 (A) 


p 2 (BK)+p z (jBK) 


7 t 2 ( ; T 


But 

So 


p 2 (BK) < p z (BK) +p z (jBK) < \\BK 


tt2 


tt2 ( : t 


td< — 

CO 


„J m) ll^ll +P(jBK) V \\BK\\ 2 -p 2 (A) 


\\BK\\ 


< — cos 

CO 


B (A) + ii (JBK) 


(61) 


(62) 


(63) 


The imaginary parts of the system poles are bounded from above by 
co = lmXj (-jA + jBKe^ J(0td ) <p (-jA) + p (jBKe~ jw ‘ 1 ) < p(-jA) +p(jBK) |cos cot d \ +p(BK) |sintau| (64) 
which can be expressed as 


® < p (—jA) + \J\ p 2 (BK) + p 2 (jBK) < p (—jA) + 11^1 


(65) 


Without loss of generality, we replace the inequality sign with equality sign for the purpose of estimation of the 
time delay margin. 

Example: Given 


A = 


0 1 

-1 1 


BK = 


0 0 
0 2 


co and t d are computed as follows: 


P(A) = 1 
P (jBK) = 0 
\\BK\\=2 

co = p (—jA) + ll-BA'H = 3 rad/sec 


1 

t d = — cos 
CO 


\\BK\\ 


= — = 0.349 sec 


Comparing the result with the exact time delay margin, the matrix measure method produces the least conser- 
vative estimation of the time delay margin. Moreover, the solution is much simpler to compute than both the Pade 
approximation and Lyapunov-Krasovskii method. 

The matrix measure method can also estimate the effective phase margin of a MIMO system. For this example, 
the phase margin estimate is 

71 

<j) = COt d = - 


4 Time Delay Margin Estimation of Adaptive Control 

4.1 Standard MRAC 

Using the bounded linear stability analysis method to approximate the standard MRAC as a bounded linear approxi- 
mation within a given time window, the time delay margin of adaptive control can be estimated by any of the methods 
previously presented. 
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1. Pade approximation: 

The first-order Pade approximation of the bounded linear approximation of the error equations and adaptive law 
are expressed as 


sE (s) = (A + B0* T <1>') E (a) - BK X E (s) + BU ad (a) - B&* r - B®* r <i>', [X m (s) - x«] 

2 — t d s 
2 + t d s 


BKr 


Xm (^) X m (^) 


2 + t d s 
2 -t d s 


2 + t d s 


+ BK, 


R(s)-R(s) 


Uad {s) = - yB T P ^ + ®J 4>; [X m (s) - E ( 5 )] 

Then the time delay margin can be computed from the following characteristic equation: 

JBB P \ 2 - jCO,t di 


det 


mi- (. 


- (A + B0* T <1>' 


BK X + B®J<Z> 


ja>i J 2 + ja>jt dj 


= 0 


( 66 ) 

(67) 

( 68 ) 


Another approach is to express the error equation in another form by differentiating the error equation and 
substituting in the adaptive law which yield 


e{t) = 


e(t) = (A + B0* t T> 


>’)e(t)- [BK X + B&J O') e (t - t d ) - yBB 1 Pe (t - t d ) + B®J o'x m (t-t d ) 

- B®* t o'i,„ ( t ) - BK X [x m (, t ) - x m ( t -t d )]+ BK r [r ( t ) - r (t - t d )\ (69) 


for t G [ti — To,tj), where = 0, tj = f,_i + Tq and i = 1,2, . . . ,n — > <». 
Then the error equation can be recast as 


z (t) = Qz ( t ) - Djz (t-t d ) + Ei (t) 


where z{t) = e(t) <?(/) and 


C,= 


A = 


E (t) = 


0 / 

0 A + B0* t O; 

0 0 
yBB T P BK X + B®J O, 

0 


[ 50, 0 ( i m (/ - t d ) - B0* O f i m (f) - //A [x m (f ) - x m (f - /</)] + BK, [r (f) -r(t-t d )\ J 


(70) 

(71) 

(72) 

(73) 


Then using the result in the previous section, the time delay margin can be found by solving the following 
characteristic equation 

JO), I -I 


det 


-^0 Q-D, ) jcOi+^I-Ci-D, 


= 0 


(74) 


2. Lyapunov-Krasovskii method: 

Applying the result of the Lyapunov-Krasovskii method from the previous section, the time delay margin can 
be computed from the following LMI: 


(Q -Difp + P (< Q - A) + (« + j3 )t d ,PD,Dj P + ^ Cj Q + ^ Dj A < 0 


(75) 


3. Matrix measure method: 

Using the matrix measure method, the time delay margin is estimated as 


t di 


— cos 

0 ), 


At (A) + m (,/A) 

II A II 


(76) 
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COi < B (- jCi ) + II Di 


(77) 


Another approach is to derive the time delay margin from the original system matrices. The characteristic 
equation of the bounded linear approximation of the error equation and adaptive law is given by 


det 


jco,I~ (a 4 B& 1 O'- j + (BK. + BQj^- 


jyBB T P 


CO; 


Then, applying the matrix measure method, the time delay margin is estimated 


0 -mu 


as 


= 0 


(78) 


1 

t(l, < — cos 
CO i 


JI (A + fi©* T <I>^ [coy \\BK x + B&J\\+C 0 i]l(-jyBB T P)] 


cof\\BK x + B&J\\ 2 +\\yBB T P\\ 2 


co;p 


j(bk x +b@J <&;) 


\yBB J P\ 


coy \\BK x + B&] 


coy \\BK x + B&J 


-\\yBB T P\\ : 


The frequency is computed as 


yBB J P 


(79) 


CO; < p 


~j(A + Be* T <p[) ] + JI [j (j BK X + B&J ) 


BK X + B&; e 


A ) o-JWdi 


— B 
co, 


co < 


B 


-j (a+B@* t 4>'^ 


<B 

|| BK X + B&J 


j(a+b &* t $;) 


(jBB J Pe~ jC0l, A 
BK X + B&J 


1 

co, ; 


yBB J P 


(80) 


1 + . 1 + 


4 1| yBB T P 


\ {Arh^A+Be^^j + llB^+BQTll} 1 


(81) 

It is noted that ft), — > as y — > Consequently, the time delay margin tends to zero as F — > This is consistent 

with the behavior of the standard MRAC. 


4.2 Scalar MRAC 

Consider an input delay scalar MRAC system with linear structured uncertainty 


x(t) = ax(t ) +b [u {t — td) + 0*x(f)] 


The reference model is given by 


The controller is given by 


x„, (0 = a m x m ( t ) + b,„r ( t ) 

u (t) = —k x x (t) + k,-r (t) — 6 (t) x (t) 
9 (f) = — rx(f) pbe(t) 


1 . The time delay margin is estimated from the matrices C, and Z), 


where 


C,= 


0 1 
0 a + b9* 


Di — 


0 0 
yb 2 p bkx + bOi 


7 = 


r 

% 



{x)dx 


(82) 

(83) 

(84) 

(85) 

( 86 ) 

(87) 

( 88 ) 
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Using the matrix measure method, the following parameters are computed analytically as 


ft), and td i are then estimated as 


a + bO* + \ (a + bd*) 2 + 1 

nQ) = ^ 

(89) 

n-jCi) = l 

(90) 

WD,)=*l P 

(91) 

II A II = \J (‘ bk x + bdi) 2 + y 2 b 4 p 2 

(92) 

~P (- jQ ) + || A|| = ^ + \j ( bk x + bQi) 2 + y 2 b 4 p 2 

(93) 


, _ 1 „„„-iM(A) + MC/A) 
"'"ft* I! All 


1 + 2-w (bk x + bOi) + y 2 b 4 p 2 


+ bO* + J(a + bO*) z + 1 + yb~p 


2J ( bk x + bQi)~ + y 2 b 4 p 2 


2. Using the original system parameters, the time delay margin is estimated as 

\bk x + bdi\ 1 


ft); = 


2 + 2 V ( bk * + bQ i) + 4 V b2 P 


td = — cos 
' ftf 


(Of ( a + bO*) \bk x + b6j\ + jb~p\J (Of ( bk x + bOj)~ + y 2 b 4 p 2 
cof ( bk x + bdi) 2 + y 2 b 4 p 2 


(94) 


(95) 


(96) 


Both approaches yield somewhat different results. The “exact” values of ft), and for the locally bounded linear 
approximation of the error equation can be determined as follows: 


det 


p j^d 

jco — (a + bO *) + (bk x + bOA e^ ]<0 ' d + ylrp 

./ft) 


= 0 


This results in two equations 

— (Of + ( bk x + b6j) ft),- sin ft),/,/, + yb 2 p cos ft),/,/, = 0 

— (a + bd* ) ft), + ( bk x + bOj) CO; cos ft),/,/. — yb 2 p sin ft),/,/. = 0 
The frequency equation is obtained as 


ft), 


(a + bO*) 2 - (bk x + bQi) 2 } co 2 - yV/? 2 = 0 


(97) 

(98) 

(99) 

( 100 ) 


The “exact” solution gives 


ft), 


\ 


(fcfcjc + fcO,-) — + + -\ / + — (bk x + bQi) + 4y 2 /? 4 /? 2 


td = — cos 

1 CQj 


cof ( a + bO *) (bk x + £0;) + (ofyb 2 p 
(of ( bk x + b6i) 2 + y 2 b 4 p 2 


( 101 ) 


(102) 
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It is well-known that t d — > 0 as as T — » °° for the standard MRAC. This behavior is exhibited in the time delay 
margin estimation by the matrix measure method and the “exact” solution since 

lim — = lim — . = 0 (103) 

00 y-^°° yj yb 2 p 


For y=0, the system is non-adaptive and the time delay margin estimation by the matrix measure method using 
the C, and Z), yields 


td = 


1 + 2 bk x 

The exact time delay margin is computed to be 

td = 


1 


+be* + ^(a+be*)~ + i 
2bk x 


a + bO * 


\J ( bk x ) 2 - ( a + bd *) 2 


bk x 


(104) 


(105) 


Example: Given a = 1, b = 1, 9* = 0.1, a m = —1, b m = 1, p = 1, 9 (0) = 0, r(t) = sin(f). The control gains are 
computed to be k x = 2 and k r - 1 . The adaptive gain is selected as T 1 . 

For the non-adaptive LTI system for which 0 (t) = 0 for all t, the time delay margin estimates and the exact value 
are shown in Table 1 . 


Method 

t d , sec 

Pade Approximation 

0.646 

Fyapunov-Krasovskii, a = p = 1 

0.137 

Matrix Measure in Part 1 

0.347 

Matrix Measure in Part 2 

0.494 

Exact Value 

0.592 


Table 1 - Time Delay Margin Estimation of Non- Adaptive System 


Thus, it can be seen that the time delay margin computed by the matrix measure method is the least conservative 
lower bound estimate of the true time delay margin among the present approaches. The difference in both approaches 
using the matrix measure method is noted. The approach using the matrices C, and l), are more conservative. 

Figure 2 is a plot of the variation of the local time delay margin estimates computed by the matrix measure method 
using the matrices C, and I), within three different time windows with 7 q = 1 sec, 7 q = 5 sec, and 7 q - 1 0 sec based 
on the bounded linear stability analysis method. It is noted that as the window size increases, the variation in the local 
time delay margin decreases. The longer time window allows a more uniform average value of the parameter y to be 
computed, thus reducing the local variation of time delay margin from one time window to the next. It appears that 
the mean value of the computed local time delay margins is relatively insensitive to the window size. In fact, the mean 
estimate of the time delay margin for the three time window sizes Tq = 1 sec, Tq = 5 sec, and Tq = 10 sec are 0.434 
sec, 0.432 sec, and 0.432 sec, respectively. 

Figure 3 is a plot of the mean value of the time delay margin estimates by all the different methods as a function 
of the unknown parameter — 1 < 9* < 1 for Tq 1 sec. Generally, 9* is not known, so in a verification setting, 
the time delay margin should be computed over all possible parameter variations within their physical bounds. Also 
plotted is the numerical evidence of the time delay margin from simulations. Comparing to the numerical evidence, the 
mean value of the time delay margin estimates by the matrix measure method using the C, and I), matrices is the most 
reasonable and non-conservative. The difference between the numerical evidence and the estimate is relatively uniform 
by about 20% to 30%. On the other hand, the second matrix measure approach over-estimates the time delay margin 
for 9* < 0 but produces a very good agreement with the numerical evidence for 9* > 0. The Pade approximation 
under-estimates for a small range of 9* between -1 and -0.5. Otherwise, it uniformly over-estimates the time delay 
margin by the largest margin among the three methods. The Fyapunov-Krasovskii method, as expected, yields the 
most conservative estimates of the time delay margin. The difference is significant. For example, for 9* = —1, it 
estimates a value of t d = 0.0356 sec as compared to the numerical value of t d = 0.719 sec which results in a difference 
of 20 times. At the other end, for 9* = —1, the Fyapunov-Krasovskii estimates t d = 0.0051 sec as compared to the 
numerical evidence of t d = 0.250 sec which is 50 times larger. Thus, in summary, the matrix measure method appears 
to give a better prediction of the time delay margin among all the three methods. 
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Fig. 2 - Time Delay Margin Estimates 



Fig. 3 - Time Delay Margin Variation with Unknown Parameter 0* 



Fig. 4 - Time Delay Margin Variation with Adaptive Gain 

Figure 4 is a plot of the mean value of the time delay margin estimates by the matrix measure method as a function 
of the adaptive gain 1 < F < 100 for To = 1 sec. It can be seen that as the adaptive gain F increases, the time delay 
margin of the adaptive system decreases as expected for the standard MRAC. 
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4.3 ^-Modification 

The (7-modification adaptive law [9] is given by 

©(f) = -r<t>(x(f)) [e T (f)Pfi + c7©(f)] (106) 

Using the bounded linear stability analysis method, the adaptive law is approximated as 

Uad ( t ) = —jB J Pe (t) - Tauad (t) + ©, <t>, [x m (f) - e (f)] (107) 

where T = A m ,„ (r), for f £ [tj — 7o, f,), where to = 0, f, = f,-_ | + To and i = 1,2, 

Then the error equation is obtained as 


e (t) = (a + - Ter) e ( t ) - (bK x + B&J <t>') e(t- t d ) + Ta (a + fi©* T <l>') e (t) 

- (yBB r P + VcjBkJ'j e(t — t d ) - TgB&* t <t>, + (t) - - ToBK x [x m (f ) - x m (t-t d )] 

+ ToBK, [ r ( t ) - r(t - t d )\ + B&J <t >)x m ( t - t d ) - B &* 1 (f) - BK X [x m (t) - x m ( t - t d )\ 

+ BK r [r(t)-r(t-t d )\ (108) 

The time margin for the (7-modification adaptive law can be estimated by the following methods: 

1 . Pade Approximation: 

The time delay margin can be found from the following characteristic equation: 


ml - (a + b©* t 4>;) 


BK X + B&J% m 

jcOi + Fa 


2. Lyapunov-Krasovskii method: 

The time delay margin for the cr -modification adaptive law can be estimated using Eq. (75) with the matrices C, 
and Dj given by 

[ 0 / 1 
Ci= To (a+B&* t <p'^) A + fi©* T <l/-rc7 <110) 


jBB 1 P + ToBK x BK x + B&J(P i 


3. Matrix measure method: 

The characteristic equation is 


det [ ml - (A + B&* 1 <) + ( BK X + B&J <lj . Ja j + JBB ' /J 
V 7 V jcOj + Ta jcOj + To 


e~ i<0i ' d i =0 (112) 


which can also be expressed as 


det /ft),/ - (a 4- B&* ' (fi'j + BK X 


B&J + TayBB 1 P + jcoj (voB&J <!>' - yBB 1 P 


cof + F 2 a 2 


o-iWdt = o 


The same approach can be applied to find the expressions for the time delay margin and frequency by working 
with the following inequalities 


0<jU A + S0* '<&,. + jU - BK, 


B&l jco i (r_oB&J<i>' l ~YBB T P 


(of+Tja 2 


co? + r 2 a 2 
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C0j < n 


-j(a + S0* t o')] +li(jBK, 

^/T ayBB Pe 


-jWd: 


F 2 g 2 




-j^ih 


®T_ 

® 2 +r 2 

ox 


co 2 + r 2 a 2 L‘ 


jj. ( yBB 1 P - FaB&J <&'■) e~ jC0i ‘ d ‘ 


(115) 


The frequency equation, however, results in a cubic polynomial which does not have an analytical solution. 

Consider the case when T — > °°. The time delay margin for the standard MRAC is known to be zero but remains 
finite for the cr-modification adaptive law. To see this, consider the characteristic equation in the limit as F — > °° 


det 


mi- (a+b©* t $;) 



yBB r P\ 

~T^~ ) 


g-mtdt 


= 0 


If r = cl where c > 0 is a constant, then 


r> 


Jp / 4> T (x(T))4>(jr(T))d/T 

aT 0 Jti-T 0 


(116) 


(117) 


is finite. 

Using the matrix measure method, the time delay margin for the cr-modification as F 
as 


— j ("a +Z?©* T <J>;') 

+ 

yBB 1 ? 
BK X + ' ^ 

L V V J 


rcr 


tdj 


M 1 

(a+r©* t <i>^ 

1 +M 1 


BK X + 7ll y/ 

) 


BK x + l 



°° can then be estimated 
(118) 

(119) 


It is clear that the time delay margin for the cr-modification remains finite even as T — > °°. It is noted that the 
frequency and time delay margin are dependent on time windows since y varies with different time windows. 


4.4 Optimal Control Modification 

The lack of robustness to unmodeled dynamics [13] and input time delay of the standard MRAC is well-known, 
as illustrated in the previous section. Increasing the adaptive gain allows the adaptation to attain a better tracking 
performance, but usually at the expense of the ability to maintain stability in the presence of unmodeled dynamics or 
time delay. To improve robustness, the two well-known robust modification methods; namely, the cr-modification [9] 
and e-modification [14], have been used extensively in adaptive control. The optimal control modification method is 
a recently developed adaptive law to address robustness due to fast adaptation [10]. This adaptive law is given by 


©(f) 


-r 


$(x(f))e T (t)PB 


vT> (x (f )) 4 > t (x (0) © (0 B J PA m l B 


( 120 ) 


where v > 0 is a weighting constant. 

In can be shown that in the limiting case as T — > °° and for a linear matched uncertainty where <f> (x (f)) = x (f), the 
optimal control modification exhibits an asymptotic linear behavior [16, 15], This behavior can be explained using the 
analysis techniques developed herein. 

The bounded linear stability analysis provides a method for analyzing an adaptive system with input time delay 
using a time windowing approach. Using the bounded linear approximation of an adaptive system, the matrix measure 
method can provide a non-conservative lower bound local estimate of the time delay margin. Unfortunately, the local 
estimation is non-unique and is dependent on the length of a time window, as is evident in Figure 2. The optimal 
control modification adaptive law has a unique property that enables it to be analyzed for stability without the use of 
the windowing approach of the bounded linear stability analysis as T — » °° for <t>(x(f)) = x(f). Thus, the time delay 
margin can be uniquely estimated for this adaptive law. 

The bounded linear approximation of the optimal control modification adaptive law is 

u a d (t) = -"/B J Pe (f) + yvB r A m J PBu ad (f) + ©^ T> ( - [x m (/) - e (f)] (121) 
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for t € [tj — To Pi), where to = 0, f,- = f,-i + To and i = 1,2, . . . ,n — > °°. 

Let G = —BB t A~ t P > 0, then the error equation is obtained as 

e (t) = (a + B0* T < $'■ - yvG^j e(t)~ (bK x + B&J <&'■) e (t - t d ) + yvG (a + B©* T <f>;) e (t) 

- (yBB r P + yvGBK^J e(t-t d ) - yvGB&* J 4>,- + 4>'x m (f ) - o'x, - yvGBK x [x m (t) - x m (t - t d )} 

+ yvGBK r [r(t)-r(t-t d )\+ B&J <I>-x m (/ - - B0* T «I>-x m (t ) - BX; V [x m (t) - x m ( t - t d )\ 

+ BK r [r(t)-r(t-t d )} (122) 

The time margin for the optimal control modification adaptive law can be estimated by the following methods: 

1 . Pade Approximation: 

The time delay margin can be found from the following characteristic equation: 

det | /ft),'/ — (a+B0* t 4>;) + [BK x +{jco,I+yvGy l + yBB 1 p ) ] (123) 


2. Lyapunov-Krasovskii method: 

The time delay margin for the cr -modification adaptive law can be estimated using Eq. (75) with the matrices C, 
and Dj given by 

r o / i 


yvG (/ 


A + B0* t 4>, A + B&* T <t>j — yvG 


yBB r P + yvGBK x BK X + B&J< I>. 


3. Matrix measure method: 


The characteristic equation for determining the time delay margin is 

detjjW- (A + B©* T <I>-) + BK x +(j(OjI+yvGy l ^j(o i B&J^' i + yBB T P^j | = 0 (126) 

Since (jcOjl + yvG) 1 is a matrix inverse, it is difficult to apply the matrix measure operator. To handle this, 
(jcOjl + yvG )- 1 can be replaced by (yft),7+ yvG)~ where G= O m j n {G) and a ndn is the minimum singular 
value. Then the characteristic equation becomes 


joy I - (a -I «©■ 1 Ay 


B&J&jCoJ + yvGyBB r P + jaij (yvGB&- T>. - yBB r P 
cof + y 2 v 2 G 2 


The solution approach is similar to that for the cr-modification. The frequency ft), does not have an analytical 
expression. 

From the fact that t d — > 0 as T — > °° for the standard MRAC, so a lower bound estimate of the time delay 
margin for which an adaptive law is stable can be estimated by the value of t d that corresponds to T — > °° or 
equivalently y — > <=°. Consider the case when y — > ° « and the matched uncertainty is linear with <f> (x (t ) ) = x(t). 
The asymptotic solution of the bounded linear approximation of the optimal control modification is obtained by 
taking the limit as y — > °° 

Bu ad (t)= l -p- l Aj n Pe(t) (128) 

Then the asymptotic error equation as y — > °° becomes 

6 0 1 ) = (a + B0* t ) e (0 - (bK x - X -p- l AlP ^ e{t-t d )~ B&* T x m (t) - BK X [x m (t) - x m ( t - t d )} 

+ BK r [r(t)-r(t-t d )) (129) 
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which is a LTI input delay equation independent of any time windowing parameters such as y. 

The characteristic equation of the asymptotic error equation can be obtained by letting y — > °° which yields 


det 


jcol - (a + B©* T ) + (bK x +^G~ 1 


-G 1 BB t P ) e 


U: 


= 0 


(130) 


where G l BB r P = -P^A^P. 

The time delay margin and frequency of the optimal control modification adaptive law as T — > °° can then be 
estimated by the matrix measure method as 


td= ~ 


JL(-jA-jB&* T)+ BK x - l -P~'AlP\\ 


At (A + BQ* J ) + AT (j [BK X - \P~ l A J m P \ ) ' 


BK X -^P-'AIP 


(-jA-jB®* J ) + 

BK X - ' p-'A J m P 

V / 

v 


The effective phase margin is estimated as 


= cos _! p(A + B®*T)+ J[ (j[BK x -lp-'Aj l P]) 
\\BK X - ±P~'AlP\ 


(131) 

(132) 


(133) 


Both the time delay and phase margins are independent of the time windows in the limit as T — » °° so the 
subscript i is dropped from the expressions. Because the time delay margin is a minimum as T — > the time 

delay margin estimate for the asymptotic solution of the optimal control modification adaptive law establishes 
a lower bound of the time delay margin for any adaptive gain T < °°. Thus, to maintain stability of the input 
delay adaptive system using the optimal control modification adaptive law, a suitable selection of the weighting 
constant v can be chosen to satisfy the time delay margin requirement and or phase margin requirement in an 
adaptive control design. In order to compute this estimate, the knowledge of the unknown parametric uncertainty 
©* must be available. 


5 Simulations 

Consider a longitudinal pitch dynamic model of an aircraft 


mV - 
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2V 
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2V 
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c, 
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0(0 

q(t) 


(134) 


where t,j is a time delay. 

A numerical model for a full-scale generic transport model (GTM) at Mach 0.8 and 30,000 ft with the flight path 
angle y = 0 is given by 


«(0 


' -0.7018 0 0.9761 


«(0 

0(0 

= 

0 0 1 


0(0 

9(0 

s. 

-2.6923 0 -0.7322 


. 9(0 
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' -0.0573 



+ 

0 




-3.5352 


$e{t-td)+ [ 0.1393 0 -0.2071 


«(0 

0(f) 
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A desired reference model of the pitch attitude is given by 


G m (t) +2£co n d m (t) + co,;d„, (t) = co,;r(t) 


(135) 


where £ = 0.85 and a>„ = 1.5 rad/sec are chosen to give a desired handling characteristic. 

Letx(f) = [ a(t) 0(f) q(t) ] T , u(t) = 8 e (t), and 0* T = [ 0« 0 0 * ] = [ 0.4 0 —0.3071 ] which 

represents a parametric uncertainty that is equivalent to an 80% reduction in the pitch damping coefficient C, llq . A 
nominal controller is designed as u nom (t) = —K x x (f) + k r r (?) where K x = y [ a$\ ( 0 „ 2£<a„ + a 33 ] = 

[ 0.7616 —0.6365 —0.5142 ] and k r = = 0.6365. The closed-loop eigenvalues are —0.6582 and— 1.2750± 

0.7902/. The nominal closed-loop plant is then chosen to be the reference model as 


’ Om(t) ' 


' -0.6582 

-0.0365 

0.9466 


’ a,„(t) ' 


' 0.0365 ' 

dm ( t ) 

= 

0 

0 

1 


dm { t ) 

+ 

0 

(f ) 


0 

-2.2500 

-2.5500 


Qm (0 

y v. 

2.2500 

Xm 



A m 


x m 


B m 


r(f) 


The control input is given by 

u (/) = — K x x(t ) +k, r (f) — © T (f)x(f) (136) 

Figures 5 and 6 are plots of the estimates of phase and time delay margins of the optimal control modification 
adaptive law for F — > °° computed by the matrix measure method from Eqs. (133) and (131) as a function of v with 
and without the parametric uncertainty 0*. Note that the phase margin generally decreases as v increases and reaches 
a steady state value, while the time delay margin reaches a maximum at about v = 1 . Thus, for practical design 
purposes, v should be kept between 0 and 1 . A large value of v produces a better time delay margin, but also results 
in a poorer tracking. For the specified uncertainty ©*, the maximum time delay margin is estimated to be 78 msec. 
Therefore, the input delay adaptive system will be stable with the optimal control modification adaptive law for any 
tj < 78 msec. 

Suppose the input time delay of the system is td = 50 msec. For this input time delay, the optimal control modifi- 
cation adaptive law produces a stable adaptation for v > 0.244. The controller is then implemented with v = 0.25 and 
r = 3000 1 selected. Figures 7 and 8 illustrate the pitch angle and pitch rate responses due to the standard MRAC and 
optimal control modification adaptive law. The MRAC with the adaptive gain F = 3000 1 does not track the reference 
pitch angle very well. High frequency oscillations are discernible in the pitch rate response. On the other hand, the 
optimal control modification adaptive law produces good tracking of the reference pitch angle and pitch rate. The 
high frequency oscillations in the pitch rate response with the standard MRAC is substantially reduced by the optimal 
control modification adaptive law. 

The time delay margin of the closed-loop system is estimated by numerical evidence to verify the lower bound 
estimate of the time delay margin for the optimal control modification. The results are shown in Table 1 . 



Time Delay Margin, msec 

No Adaptation 

550 

MRAC 

50 

OCM 

100 


Table 1 - Time Delay Margin Estimates 

The numerical evidence of the time delay margin for the optimal control modification adaptive law is estimated to 
be 100 msec. This is a factor of two larger than the time delay margin of 51 msec as estimated by the matrix measure 
method for v = 0.25. 
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Fig. 5 - Phase Margin Estimation of Optimal Control Modification 



Fig. 6 - Time Delay Margin Estimation of Optimal Control Modification 



Fig. 7 - Pitch Angle Response 
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Fig. 8 - Pitch Rate Response 


6 Conclusions 

A method for analyzing input delay adaptive systems is presented. The bounded linear stability analysis approximates a 
nonlinear adaptive system as a bounded linear approximation within a time window. The windowing analysis is used to 
analyze local stability of the bounded linear approximation to estimate local stability behavior of the original nonlinear 
system within a given time window. Analytical approaches for computing the time delay margin are presented for three 
different methods: Pade approximation, Lyapunov-Krasovskii method, and matrix measure method. The windowing 
analysis results in non-unique estimates of the time delay margin of the bounded linear system since the results are 
dependent on the time window size. The non-uniqueness is a drawback of the windowing analysis. It is found, 
however, the mean value of the local time delay margins are relatively insensitive to the time window size. Among the 
three different methods, the Pade approximation using a first-order Pade rational polynomial generally tends to over- 
predict the time delay margin, while the Lyapunov-Krasovskii tends to be highly conservative in the time delay margin 
estimation. Neither of these estimations may be deemed practical for adaptive control design purposes. The matrix 
measure method seems to be able to predict the time delay margin with a reasonable agreement with the numerical 
evidence. Moreover, the method is much simpler to use and does not require solving a linear matrix inequality as in 
the case of the Pade approximation or Lyapunov-Krasovskii method. 

The asymptotic behavior of the time delay margin as the adaptive gain tends to infinity is studied. The standard 
model-reference adaptive control has zero time delay margin as the adaptive gain tends to infinity, as expected. The 
time delay margins for both the cr-modification adaptive law and optimal control modification adaptive law remain 
finite as the adaptive gain tends to infinity. The optimal control modification adaptive law also exhibits another useful 
feature in that the asymptotic value of the time delay margin is independent of the time window and the closed-loop 
input delay adaptive system tends to a LTI system. This behavior enables a lower bound of the time delay margin to be 
estimated with ease using the matrix measure method to guarantee stability for the input delay adaptive system. Flight 
control simulations demonstrate that the time delay margin estimation by the matrix measure method provides a good 
lower bound. 
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